Relaxation and Pre-thermalization in an Isolated Quantum System 
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Understanding relaxation processes is an important unsolved problem in many areas of physics. 
A key challenge in studying such non-equilibrium dynamics is the scarcity of experimental tools for 
characterizing their complex transient states. We employ measurements of full quantum mechanical 
probability distributions of matter-wave interference to study the relaxation dynamics of a coher- 
ently split one-dimensional Bose gas and obtain unprecedented information about the dynamical 
states of the system. Following an initial rapid evolution, the full distributions reveal the approach 
towards a thermal-like steady state characterized by an effective temperature that is independent 
from the initial equilibrium temperature of the system before the splitting process. We conjec- 
ture that this state can be described through a generalized Gibbs ensemble and associate it with 
pre-thermalization. 



Despite its fundamental importance, a general under- 
standing of how isolated quantum many-body systems 
approach (thermal) equilibrium is still elusive. While 
theoretical concepts such as the quantum ergodic theory 
or the eigenstate thermalization hypothesis [IHSj try to 
infer requirements for a system to be able to undergo 
relaxation, it is still unclear on what timescale this will 
happen. Most prominently, in situations where conser- 
vation laws or the presence of many constants of mo- 
tion inhibit efficient relaxation, many-body systems are 
expected to display a complex behavior. An intriguing 
phenomenon which has been suggested in this context 
is pre-thermalization 4|. It predicts the rapid establish- 
ment of a quasi-stationary state that differs from the real 
thermal equilibrium of the system. Full thermalization, 
i.e. relaxation to the real thermal equilibrium, if present 
at all, is expected to occur on a much longer time scale. 
Pre-thermalized states have been predicted for a large va- 
riety of physical systems [SHE] and it is conjectured that 
they can be described by equilibrium statistical mechan- 
ics through a generalized Gibbs ensemble [U |3j [9] . Here 
we present a direct observation of such a state. 

Systems of ultracold atoms provide unique opportuni- 
ties to experimentally study such non-equilibrium prob- 
lems because of their almost perfect isolation from the 
environment. Moreover the timescales for internal re- 
laxation processes (collisions) are easily accessible in ex- 
periments. Consequently, there have recently been vari- 
ous studies about non-equilibrium dynamics in ultracold 
atom systems [TTH - HS] . 

One-dimensional (Id) Bose gases are of particular in- 
terest because they inherently contain strong fluctuations 
and dynamics: At finite temperature many longitudinal 
modes of the system are populated which manifests itself 
in the rich spatial structure and dynamics in their local 



phase. This is in stark contrast with three-dimensional 
condensates, where the existence of long-range order al- 
lows to characterize the state with a single, global phase. 
In addition a homogeneous Id Bose gas is a prime ex- 
ample of an integrable quantum system |17j . The near- 
integrability of experimentally realized trapped Id Bose 
gases thus opens up the possibility of studying dynamics 
and relaxation close to an integrable point. 

In our experiment (Fig. 1) we start from a single Id 
Bose gas of 87 Rb in the quasi-condensate regime [TH] pre- 
pared in an elongated microtrap on an atom chip |19j . We 
prepare the initial state for our evolution by rapidly and 
coherently splitting the single Id gas, producing a system 
of two uncoupled Id Bose gases in a double-well poten- 
tial. The two trapped Id gases only differ by the quantum 
shot-noise introduced in the splitting (Figs. la). They 
have almost identical longitudinal phase profile, and are 
therefore strongly correlated in their phase. In con- 
trast, two independently created quasi-condensates have 
vastly different and uncorrelated phase profiles (Fig. 1 
right column) . The strongly correlated phase of the two 
gases after splitting reflects the memory that they orig- 
inally come from a single quasi-condensate. Our experi- 
ment studies how this memory about the the initial state 
evolves, decays in time, and if a thermal equilibrium state 
corresponding to two independent and classically sepa- 
rated quasi-condensates is reached in the evolution. 

We probe the evolution of the local phase difference 
between the two quasi-condensates by matter-wave in- 
terference (Fig. lc) . The system is allowed to evolve in 
the double-well for some time t e before the two Id gases 
are released from the trap and allowed to interfere in 
time-of-flight. The local phase difference along the axial 
length of the system directly translates to a shift of the 
interference peaks (Fig. lc) . To probe the strength of the 
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FIG. 1: Experimental scheme, (a) An initial phase fluctuating Id Bose gas is split into two uncoupled gases with almost 
identical phase distributions t/>i(«) and ^2(2) (represented by the black solid lines) and allowed to evolve for a time t e . (b) 
At t e = ms, fluctuations in the local phase difference A<f)(z) between the two gases are very small and the corresponding 
phase correlation length is very large. During the evolution these relative phase fluctuations increase and the correlation length 
decreases. The main question we address in this paper is, if or when this system will reach the corresponding thermal equilibrium 
of uncorrelated phases as characterized by the initial temperature T, and thermal coherence length At. In experiment, this 
situation can be prepared on purpose by splitting a thermal gas and cooling it into two independent gases |16| . (c) shows 
typical experimental matter-wave interference patterns obtained by overlapping the two gases in time-of-flight after different 
evolution times. Differences in the local relative phase lead to a locally displaced interference pattern. Integrated over a length 
L, the contrast C(L) in these interference patterns is a direct measure of the strength of the relative phase fluctuations, (d) 
Due to the stochastic nature of the phase distributions, repeated experimental runs yield a characteristic distribution P(C 2 ) of 
contrasts, which allows one to distinguish between the initial state, an intermediate pre-thermalized state and the true thermal 
equilibrium of the system. 



fluctuations in the local phase difference A(j)(z), we inte- 
grate the interference pattern longitudinally over a vari- 
able length L and extract from the resultant line profile 
our main experimental observable: the integrated con- 
trast C(L) (Fig.lc). For the initial state, the local phase 
difference is close to zero everywhere along the quasi- 
condcnsates, and thus the integrated interference con- 
trast C(L) is large for all integration lengths L. During 
the course of the evolution, the phase difference varies in 
the longitudinal direction due to the strong fluctuations 
inherent in Id systems, which results in the reduction of 
C(L) starting with long integration lengths. Thus the 
measurements of C(L) allow the characterization of the 
unique dynamics of Id quasi-condensates. 



The mean squared contrast (C(L) 2 ), similarly to the 
coherence factor used in [TT] , is a direct measure of the 
integrated two-point correlation function of the relative 
phase between the two halves of the system [20j [2Tj . Inte- 
grating over the whole length of the interference pattern 
we observe (Fig. 2a) an initial rapid decay in (C 2 ) on a 
time scale of approximately 10 ms, after which a quasi- 
steady state emerges which slowly evolves further on a 
second, much slower time scale. 

The initial rapid decay in Fig. 2a is analogous to the 
one observed in the experiment presented in |llj which 
was limited to evolution times t e < 12 ms by longitu- 
dinal dynamics introduced in the splitting. Significant 
improvements in the experimental techniques (see SOM) 
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FIG. 2: (a) Evolution of the mean squared-contrast (C 2 ) for 
interference patterns integrated over the whole length of the 
Id systems. We observe a rapid decay followed by a long 
slow further evolution. Error bars are standard errors of the 
mean, inset: Experimental non- equilibrium distributions of 
C 2 1 (C 2 ) at t e = 22 ms, 47 ms, and 122 ms respectively (his- 
tograms) and a fit of a theoretical equilibrium distributions 
leading to T e g = 15±| nK, 14±| nK, and 31±g nK respec- 
tively (red solid line). For comparison the calculated equilib- 
rium distributions for T — 78 ± 10 nK (blue dashed line) are 
added, (b): Evolution of T e ff for the whole data set extracted 
by fitting equilibrium distributions. A linear fit indicates an 
increase of T e g over time of 0.14 ± 0.04nK/ms. The yellow 
area indicates the measured heating rate of our atom trap of 
0.11 ± 0.06 nK/ms 



allowed us now to reveal the long time behavior. In the 
following we will first show that this steady state is not 
the expected thermal equilibrium and associate it with 
pre-thermalization [J. 

To probe the nature of this quasi-steady state we start 
by employing tools developed to characterize equilibrium 
systems [20 23 and capture higher-order correlations in 
the system through the higher moments (C 2n ). For 
this we extract the full quantum mechanical probability 
distribution function (FDF) P(C 2 )dC 2 , which gives the 
probability to observe a value C 2 in the interval between 
C 2 and C 2 +dC 2 . The higher moments (C 2n ) are directly 
related to P(C 2 ) by (C 2n ) = J C 2n P(C 2 )dC 2 . Conse- 



quently the FDF is a direct measure of all even relative 
phase correlations between the gases and hence deter- 
mines the state of the system with unprecedented detail 
[5ni[2I]- In particular, high phase coherence between the 
two halves of the system results in a peaked Gumbel-like 
distribution, whereas the distribution is exponential in 
form when the phase coherence is low [2"0Tf2"3"] . 

Using a statistically large set of data we can map the 
time evolution of the FDFs for different length scales L. 
For times > 12 ms, i.e. directly after the initial rapid evo- 
lution shown in Fig. 2a, we find remarkable agreement of 
the measured FDFs with theoretical equilibrium distri- 
butions. We extract an effective temperature T e g from 
a simultaneous fit to the measured data on all length 
scales probed (see insets in Fig. [2^,). Surprisingly, im- 
mediately after the fast decay at t e = 12, 17, 22 ms we 
find: T c g — 13±|, 13±|, 15±| nK respectively, which 
is more then a factor of five lower than the initial tem- 
perature of the un-split system (T = 78 ± 10 nK). The 
observed steady-state hence cannot be the true thermal 
equilibrium state of the system. (For a direct comparison 
of FDFs for the thermal equilibrium distributions in the 
same double-well system, see SOM.) 

In contrast, for evolution times t e < 12 ms the shapes 
of the measured FDFs are not consistent with equilib- 
rium theory. The thermal-like appearance of the state is 
established only during the evolution of the system. 

To analyze the subsequent further slow evolution ob- 
served in Fig. [2^,, we extract the effective temperature for 
all times after the initial decay. The measured values of 
T e g are plotted in Fig.[2]D. We find an increase of T e g 
over time of 0.14±0. 04 nK/ms. This is, however, consis- 
tent with the measured heating rate of our atom trap of 
0.11 ± 0.06 nK/ms which we characterized independently 
using equilibrium quasi-condensates (see SOM). This in- 
dicates that either no thermalization is present in this 
nearly integrable system, or, if it is present, that it is a 
very slow process. 

To describe the fast evolution from the splitting to the 
emergence of the quasi-steady state, we employ a fully 
integrable theory based on a Tomonaga-Luttinger liquid 
formalism f2U (for details see SOM). The evolution 
of the local phase difference between the two halves of 
the system A<j>(z) is thereby described by a set of un- 
coupled collective modes with momentum k, i.e. sound 
waves, which modulate the relative density and phase at 
a wave-length A = 2-7r/fc and with an amplitude given 
by the population of the mode. A sudden splitting cre- 
ates an equipartition of energy between all the fc-modes, 
which initially are all in phase (SOM). The rapid evo- 
lution of the system seen over the first ~ 10 ms is then 
the dephasing of these /c-modes. The FDFs calculated 
by this integrable theory f2H HSJ , using input parameters 
independently extracted from the experiment, show re- 
markable agreement without any free parameter (Fig.|3|. 

The model also predicts a steady state, to which the 
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FIG. 3: The evolution of the squared-contrast distribution 
P(C 2 ) for different integration lengths L. Experimental data 
is plotted using histograms and the theoretical simulations 
as (red) curves. For integration over the full cloud length, 
the distribution rapidly becoming exponential with increasing 
evolution time t e . For the shortest integration length, the 
distribution preserves a non-zero peak showing a persisting 
memory of the correlations of the initial state |26) . The light 
(pink) shaded areas denote the errors of the experimentally 
measured theory input parameters. 



integrable system will relax: The dephasing, along with 
the equipartition of energy between the /c-modes intro- 
duced by the fast splitting, results in the FDFs of the 
quasi-steady state being indistinguishable from those of 
a system in thermal equilibrium at some effective tem- 
perature T e ff, which is determined by the energy given to 
the relative degrees of freedom by the quantum shot noise 
introduced in the splitting. The full calculation gives [55] 



kBT eS = gp/2, 



(1) 



where g — 2huj±a s is the Id interaction strength for 
particles with scattering length a s trapped in a tube 
with transversal confinement lj±_ and p is the Id den- 
sity of each half of the system. For the parameters 
used in the data presented in Fig. [3] the model predicts 
T e ff = 11 ± 3 nK, in very good agreement with our obser- 
vations of T cff = 14±4nK, 17±5 nK, and 14±4nK for the 
evolution times of 12 ms, 17 ms, and 27 ms, respectively. 

Moreover our integrable model predicts (Eq. [I]) that 
the effective temperature should be linearly dependent 
on the initial Id density, and independent of the initial 
temperature. Both these predictions are confirmed by 
extending the experiments over a wide range of initial 
conditions (Fig. [1]). 

The apparent systematic offset of the experimentally 
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FIG. 4: a) Dependence of T c s on p and b) Independence 
of T c ff from the initial temperature T of the system before 
splitting, corrected for the scaling of T e g with density. The 
colors encode different datasets. The (black) solid line cor- 
responds to the theoretical prediction fcsT c fi = gp/2. The 
black (green) data point in a and b correspond to the dataset 
presented in Fig. 3, (Fig. 2) respectively. 



derived T c g and the theoretical prediction in both Figs. 
[4^, and[4jD can be attributed to imperfections in the exper- 
imental splitting process, which in the model is assumed 
to be instantaneous. This finite-time splitting is also the 
reason that the agreement between the experiment and 
theory in Fig. [3] is less good for very early times. 

Nevertheless the first milliseconds of the observed dy- 
namics are well-captured by the integrable Luttinger liq- 
uid theory. The large number of conserved quantities 
in this integrable system prevents thermalization. Our 
experimental realization of a Id system is however not 
completely integrable and will eventually thermalize. 

Dynamics beyond the harmonic Luttinger Liquid 
model is required to couple symmetric and anti- 
symmetric modes and give rise to thermalization [27] . 
A mechanism that is expected to come into play is inter- 
actions of particles that go beyond two-body collisions, 
like three-body processes connected with higher radial 
trapping states [281130] . In our present experiment, these 
processes that can lead to full thermalization are much 
slower than the de-phasing of the collective modes, and 
thus allow the clear observation of the dynamics domi- 
nated by the close-by integrable system. It is of great 
interest to investigate the physics of thermalization in 
the future and study, for example, how far away from 
integrability one has to go to see full thermalization and 
probe its time-scale. 

In view of our present analysis, the observed decay 
of the coherence factor in the experiment of Hofferberth 
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et al. [TT] has to be reinterpreted. In agreement with 
our present experiment it shows the same fast 'inte- 
grable' de-phasing of relative modes in the spit Id sys- 
tem [2H HHJ [3T] , and not full decoherence and thermal- 
ization as originally interpreted by comparison with the 
theoretical description of Burkov et al. [27J (for more de- 
tails see SOM) We point out that for the present exper- 
iment, even independent of our theoretical model, the 
observed independence of T e g from the initial temper- 
ature provides direct experimental evidence that we do 
not observe thermalization. 

In summary, the quasi-steady state found in our ex- 
periments provides the first direct observation of pre- 
thermalization, as predicted to appear in non-equilibrium 
systems close to an integrable point We would like 
to point out that we observe pre-thermalization in the 
relative degrees of freedom. The initial thermal energy 
is still stored in the common mode fluctuations of the 
two halves of the system which is are not probed by the 
interference pattern. 

The effective temperature being significantly different 
from the kinetic temperature supports the prediction 
that such a state requires a description by a generalized 
Gibbs ensemble HJ HE]. 

We note that the thermal-like appearance of the pre- 
thermalized state is not necessarily generic but a special 
property of our system and is due to the fast splitting 
process which puts equal energy into all (relative) collec- 
tive modes of the system. 

Our experiment also directly shows that the two sepa- 
rated many-body systems retain memory of their initial 
state for a time much longer than the randomization of 
the global phase would suggest, and that genuine de- 
coherence which would erase the memory did not yet oc- 
cur, i.e. the two Id systems did not yet emerge as two 
classically separated entities. 

The timescale over which this pre-thermalized state 
persists remains an open question. It is directly related 
to the open problems of how two quantum-mechanically 
entangled objects reach classicality, the properties of the 
hypothetical quantum KAM theorem [I], and the very 
nature of thermalization itself. 
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